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We analyse the impact of systematic effects due to the scale dependence of QCD corrections 
in combination with the use of different jet clustering algorithms in the measurement of 
the mass in the fully hadronic decay mode of W + W~ pairs produced at LEP2. We 
consider higher order contributions induced by both virtual and real gluon radiation onto 
the electroweak CC03 and CC11 channels through 0(a s ) at the parton level. We find 
that the associated uncertainties can be of order 100 MeV, thus competitive with those 
possibly arising in the non-perturbative regime and indeed above the current experimental 
estimates. 
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1 Introduction 



Over the past years, LEP2 has been producing and studying bosons. One of the main goals 
of such a collider was the determination of Mw± with a target accuracy of 40-50 MeV. This has 
apparently been achieved. By combining their data in all possible W ± decay channels, the four 
LEP experiments quoted the following result: 

M w ± = 80.427 ± 0.046 GeV, (1) 

that compares rather favourably with the estimates obtained at pp colliders [|TJ. This measure- 
ment is extremely important: if combined with an improved determination of the top mass, 
m t (soon to be performed at the Tevatron, during Run 2), it can lead to a rather stringent 
prediction of the Higgs mass, from a fit to high precision electroweak (EW) data. 

One of the experimental strategies adopted to measure the W ± mass at LEP2 has been the 
kinematic reconstruction of the resonance through the momenta of its decay products^], 
e.g., in the fully hadronic channel: e + e~ — > W + W~ — > 4 jets. A 'cleaner' measurement is 
certainly performed in the semi-leptonic channel, i.e., e + e~ — ► W + W~ — > 2 jets I plus an 
undetected neutrino (with £ = e,fi). However, the contemporaneous presence in this case of 
missing energy in the final state and of photon radiation in the initial state (ISR), loosens the 
kinematic constraints that can be applied in the W ± mass reconstruction procedure. Besides, the 
fully hadronic decay rate has a somewhat higher statistics than the semi-leptonic one. Therefore, 
although the event reconstruction is made harder in multi-jet final states by the larger number 
of tracks in the detector and by the usual uncertainties related to measuring jet energies and 
directions (a task much less complicated in the case of leptons), the W + W~ — > 4 jets mode 
represented an accurate means of determining M\y± at LEP2 0. In fact, as shown in Ref. [@, 
the separate results obtained from the W + W~ — > qq'QQ' and W + W~ — > qq'ivi channels are 
consistent, with a difference in mass which is very small: 

AM W ± (qq'QQ' - qtftvt) = +9 ± 44 MeV. (2) 

Nonetheless, systematics errors on M w ± are somewhat larger in the qq'QQ' than in the qq'ive 
channels: see Table 2 of Ref. [0. 

Let us examine then more closely the kind of problems associated with the W + W~ — » 
4 jets signature. One of the issues is the problem of estimating theoretical biases due to the 
relatively unknown 'colour-reconnection' (CR) || and 'Bose-Einstein correlation' (BEC) |D 
effects (see || for a theoretical review and || for an experimental one). Things go as follows. In 
e + e~ — > W + W~ — ► 4 jets, one should expect some interference effects between the two hadronic 
W ± decays, simply because the decay products from the two different gauge bosons can overlap 
considerably in space-time. In fact, at LEP2 energies, the separation between the two 
decay vertices is ~ 0.1 fm, that is, much smaller than the typical hadronisation scale, ~ 1 fm. 
Hence, the two hadronic W ± decays can no longer be considered as separate, since final-state 
interactions (CR) and/or identical-particle symmetrisation (BEC) can play a non-negligible 

2 An alternative method is the so-called 'threshold scan', wherein a value for M w ± is fitted to the shape of 
the e + e~ — » W + W~ total cross section for yfs in the vicinity of 2M W ± (the result quoted in eq. (|l|) does also 
include measurements obtained in this way). 



1 



role, possibly leading to an apparent 'shift' in the reconstructed W ± mass resonance || [|. 
Unfortunately, because of our current lack of understanding of non-perturbative QCD|, such 
interference effects can only be estimated theoretically in the context of different 'models'. Most 
of the latter can be constrained by looking at experimental observables which are sensitive to 
either phenomenon. For example, CR effects would be manifest in the central region, particularly 
in the single-particle distributions at low momenta, as produced in the fully hadronic versus the 
semi-leptonic channel, whereas BECs would lead to an increase in the correlation function for 
W + W~ , as compared to that for a single ||11||. While the jury is still out in the case of CR 



effects, there is an increasing evidence that BEC effects are very small if at all present In 
practice, the values that the LEP experimental collaborations assign to the systematic errors on 
Myy± due to CR and BEC effects range from 30 to 66 MeV and from 20 to 67 MeV, respectively 
(40 and 25 MeV are adopted in the combined results) 0. 

Other problems in the fully hadronic decay channel of W + W~ pairs are associated with the 
definition of 'jets'. The problematic here is twofold. Firstly, because two identical decays take 
place in the same event, one has the phenomenon of mis-pairing of jets. That is, even in the 
ideal case in which all tracks are correctly ascribed to the parton from which they originate, 
one has to cope with the ambiguity that it is in practice impossible to uniquely assign any pair 
among the four reconstructed jets to the parent W ± on the sole basis of the event topology. Of 
all possible combinations of di-jet systems, only one is correct. Thus, an 'intrinsic' background 
exists in W + W~ — > 4 jets events, in terms of simple combinatorics. Secondly, because of the 
large hadronic multiplicity, one also has the phenomenon of mis-assignment of tracks. In this 
case, the ambiguity stems from the fact that a track assigned to a jet, the latter eventually 
identified as a parton originating from one of the W ±, s, might have actually been produced in 
the fragmentation of another parton coming from the second decay. 

Both these phenomenological aspects are clearly dependent upon the 'jet clustering algo- 
rithm' (see Ref. [|13|] for a review), wherein the number of hadronic tracks is reduced one at a 
time by combining the two most (in some sense) nearby ones (hereafter, we will quantify the 
'distance' between two particles % and j by means of a variable denoted by yij). This (binary) 
joining procedure is stopped by means of a resolution parameter, y CVLt , and the final 'clusters' 
yielding yij values all above y cut are called jetsf]. In Ref. fll3| , it was precisely this dependence 
that was investigated, by using standard Monte Carlo (MC) simulation programs based on a 
parton shower (PS) approach (see |TJ), such as HERWIG [|TJ, JETSET/PYTHIA |T§ and 



ARIADNE fl7| . The results presented there did show a rather dramatic effect in the recon- 
structed M w ± values, due to the choice of the jet finder, and of its resolution parameter as 
well. However, any shift on M w ± of this sort can be estimated accurately, as it is simply due 
to kinematic effects induced by the jet clustering algorithm itself in reconstructing the quark 
momenta starting from the PS (or after hadronisation). In practice, it can be treated as a well 
quantifiable correction to be applied to the reconstructed M w ± value, in order to reproduce 
the true one. Using the same MC programs, one can also determine the typical size of the 



3 In fact, the perturbative effects of CR are expected to be small, because of order ~ (Cfc*s) 2 /Nq x T\y± /M^i 
jjj - see also Refs. ^ |9| - so are partonic 'Fermi-Dirac correlations' (as opposed to hadronic BECs) pfj| . 

4 Here and in the following, the word 'cluster' refers to hadrons or calorimeter cells in the real experimental 
case, to partons in the theoretical perturbative calculations, and also to intermediate jets during the clustering 
procedure. 
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systematic errors due to the hadronisation process, by comparing the outputs of the various 
programs. Finally, background effects can be accounted for by exploiting the numerous event 



generators available on the market for e + e — * 4 jets, both in EW |18| and QCD [|T^]-[|T^], ]19| 
processes (see Ref. |20| for a dedicated study of the impact of such QCD background effects in 
e+e- W+W~ -f 4 jets). 

Other systematic effects remain instead quite beyond control. These are intimately related 
to the way predictions are made within standard perturbation theory. That is, to the fact 
that only a finite number of terms of a perturbative series are generally computed over all the 
available phase space. Or alternatively, that only over a restricted region of it, all terms of a 
series can be summed to all orders. Whereas the availability of the latter is in general more 
crucial to the estimation of a total cross section, that of the former can be decisive for the study 
of more exclusive observables. Given the relative size of the EW and QCD coupling 'constants' 
at LEP2, it is clear that the dominant higher order effects will be due to the emission and/or 
absorption of gluons. 

Several QCD effects entering e + e~ — > W + W~ —>■ 4 jets events have been studied so far. For 
a start, it should be mentioned that the amplitude for e + e~ — > W + W~ —>■ qq'QQ' (the so-called 
CC03 channel) is quite trivial to derive, in fact, more of a textbook example. It represents 
the lowest-order (LO) contribution to the e + e~ — > W + W~ — > 4 jets signal. Higher-order QCD 
contributions involving gluons are, for example, the real ones (i.e., tree-level processes): e + e _ — > 
W + W~ — > qq'QQ'g and e + e~ — > W + W~ — > qq'QQ'gg events, which have been calculated 
in Refs. [21J and ||, respectively, as well as e + e~ — > W + W~ — > qq'QQ'g*, with the gluon 



splitting in a quark- ant iquark pair, which was considered in Ref. [|10| (see also [p2]). One-loop 
QCD corrections to e + e~ — > W + W~ — > qq'QQ' are also known to date [F£J, and they have 
been interfered with the LO amplitudes and eventually combined with the single real gluon 
emission contribution of Ref. |2T[ into the complete 0(a s ) result In Ref. [53|, the full 



0(a s ) corrections were computed for the case of the so-called CC11 channel fT8| |, also including 



irreducible background effects in addition to e + e~ — > W + W~ production and decay. Finally, 
two-loop effects due to the virtual exchange of two gluons between the two quark pairs in 
hadronic W + W~ decays were estimated in Ref. @ in the 'soft limit' and found to be either 
small (colour- singlet exchange) or large (colour-octet exchange) but symmetric around Myy±, 
hence unobservable in general. 

We make use here of the calculations of Refs. in order to assess the size of the 

typical theoretical error due to the truncation of the perturbative series at order a s and the 
systematic effects that it introduces in observable quantities, primarily, in the 'line-shape' of 
the W+ resonance, as determined by using different jet finders to select the hadronic sample. 
Our motivations to carry out such a study are dictated by the following considerations. For a 
start, NLO corrections to both CC03 and CC11 have been found to be rather large in general 



f23| , p4| , with their actual size clearly depending upon the algorithm used. Furthermore, it is 
well known that differential distributions are typically more sensitive (particularly in presence of 
cuts over the phase space available to gluon emission) to higher order effects than fully inclusive 
quantities, such as total cross sections, where virtual and real contributions tend to cancel to a 
larger extent. Besides, in the case of the W^-mass line-shape, one would expect the distortion 
effects to be mainly induced by relatively hard and non-collinear gluons, which should be better 



modelled by an exact NLO calculation than by the PS models exploited in Ref. [13 
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2 Results 



All QCD predictions have an intrinsic dependence on an arbitrary scale, hereafter denoted by 
fi, entering at any order in a s . This scale is not fixed a priori. On the one hand, although the 
structure of the QCD perturbative expansion does not prescribe which value should be adopted 
for /i, an obvious requirement is that it should be of the order of the energy scale involved 
in the problem: i.e., the CM energy y/s (see Ref. [^5] for detailed discussions). On the other 
hand, the physical scales of gluon emissions that actually give rise to multi-jet configurations 
are to be found down to the energy scale ^y C utS- In practice, one should avoid building up large 
logarithmic terms related to the (unphysical) 'mismatch' between the process scale ~ y/s and 
the emission scale y/y cu tS and it is well known that it may be necessary to adopt a different 
scale for each observable in order to best describe experimental data taken at fixed y/s [p6| . It is 
precisely the /i-dependence of the truncated perturbative series that is treated very differently by 
each jet-clustering algorithm |H| |27|, ^HJ and the corresponding effects on observable quantities 
are what we aim to study. In our calculation, QCD effects appear through 0(a s ) only, so that the 
//-dependence is merely the one affecting the strong coupling constant at lowest order. Whereas 
its impact is trivial to assess (and account for) in the case of total inclusive rates, this is no longer 
true for differential quantities (such as mass distributions), because of the different kinematics 
of lowest and higher order contributions, respectively. Notice, however, that a problem arises 
when studying the scale dependence of a s results for algorithms based on different measures, as 
for the same y cn t the total cross section at NLO can be significantly different. A more consistent 
procedure was outlined in Ref. [P^j : that is, to compare the NLO scale dependence of the various 
schemes not at the same y cnt value, rather at the same LO rate. This is our approach. 

In order to make more manifest the effects of the interplay between the 0(a s ) corrections 
to W + W~ — > 4 partons events and the jet clustering schemes tested, we have not considered 
here Coulomb corrections to CC03 |29| (their relevance is anyway modest beyond the W + W~ 
threshold). We also have neglected the implementation of the mentioned BEC and CR phe- 
nomena, as these mainly arise in the non-perturbative domain. Similarly, hadronisation and 



detector effects were not investigated, nor those due to ISR. We refer the reader to Refs. JT3 
and HU, respectively, where their impact was studied in detail. 

As centre-of-mass (CM) energy representative of LEP2 we have used the value \/s = 175 
GeV. As for the parameters of the theory, we have adopted (in the fixed-width approach) 
M z o = 91.189 GeV, T z o = 2.497 GeV, M w ± = 80.430 GeV, T w ± = 2.087 GeV, sin 2 ^ = 0.231, 
a em = 1/128.07 and the one-loop expression for a s (for consistency), with Aq£= 4 = 0.283 GeV 
(yielding a s (y/s) = 0.123)0. Furthermore, we have kept all quarks massless as a default, in order 
to speed up the numerical evaluations. Electron and positron have mass zero too, so has the 
neutrino. Also, we have neglected Cabibbo-Kobayashi-Maskawa (CKM) mixing terms (i.e., we 
have taken the CKM rotation matrix to be diagonal). 

As jet clustering schemesf], we have used a selection of the binary ones, in which only two 



5 The one-loop values of a s adopted here are consistent with the two- loop one extracted from experimental 
fits to shape variables at the Z° peak, a s (Mz«) = 0.116, for the same choice of Aqqq 4 . 

6 We acknowledge here the well admitted abuse in referring to the various jet 'finders' both as algorithms and 
as schemes, since the last term was originally intended to identify the composition law of four-momenta when 
pairing two clusters: in our case, the so-called E-scheme, i.e., = +Pj (other choices have negligible impact 
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e+ 4 jets (CC03) 




Figure 1: NLO cross sections for e + e~ — > W + W~ — > udsc (CC03), as a function of the QCD 
scale /i, as given by four different jet clustering algorithms, at y/s = 175 GeV. (The values of 
y cnt are chosen such that the LO rates are approximately equal for all schemes.) Total hadronic 
rates are obtained by multiplying those above times four. 



objects are clustered together at any step. These are the following. The JADE (J) one [P0| , 
which uses as a measure of separation (or 'metric') the quantity 

j _ 2E i E j (l-cos9 ij ) 
Vij - - • {J) 

The Durham (D) |3]] and the Cambridge (C) |27j ones, both using 

d c 2min(^,^)(l-cos%) 
Vij = Vij = '- • (4) 

(The Cambridge algorithm in fact only modifies the clustering procedure of the Durham jet 
finder.) We also have adopted the LUCLUS or LUND (L) jet finder ||32|| , for which one has 

l = 2\Pi\ 2 \Pj\ 2 (1 -cos eg) 



on our conclusions) 
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however, with the same clustering procedure of the Cambridge scheme and without 'precluster- 
ing' and 'reassignment' (see Ref. f32||), i.e., as done in Ref. QT3]| (where it was labelled as CL). 



In eqs. (|3|)-(|5]), ^i(\Pi\) an d Ej(\pj\) are the energies(moduli of the tree- momenta) and % the 
angular separation of any pair ij of particles in the final state, to be compared against the res- 
olution parameter y cnt - The choice of these particular schemes has a simple motivation. The D, 
C and L ones are different versions of 'transverse-momentum' based algorithms, whereas the J 
one uses an 'invariant-mass' measure (the numerator of eq. (^) coincides with the invariant mass 
of the partons ij, when the latter are massless, as is the case here). In fact, these two categories 
are those that have so far been employed most in phenomenological studies of jet physics in 
electron-positron collisions, with the former gradually overshadowing the latter, thanks to their 
reduced scale dependence in higher order QCD, e.g., in the case of the 0(a 2 s ) three- JTJl, [28|, 



and 0(a^) four-jet rates [33], and to smaller hadronisation effects in the same contexts [13, 28 



Our multi-jet sample is selected at the parton level, by requiring a final state with at least 
four resolved objects (i.e., with all their y^s above a given y cn t)- When five objects survive, 
the two yielding the smallest yij value (according to the metric used) are joined together, so 
to always produce a four-particle final state. In doing so, we conform to typical experimental 



approaches: see, e.g., Ref. [3j|. The impact of a different treatment of five-jet contributions was 



assessed in Ref. [23 



Fig. |l| illustrates the dependence of the CC03 NLO rates upon the unknown /x scale, for our 
four default jet clustering algorithms, for representative choices of y cut such that the LO rates 
in the various schemes are approximately equal. (In fact, differences at LO are typically within 
15%; the actual numbers being: o LO = 1.38(1.32)[1.32]{1.54} pb for J(D)[C]{L}.) The variation 
of the NLO rates with /i, for values of the latter ranging between 5 GeV and i/iQ, denoted by 
^onlo^)? depends upon the jet algorithm, varying significantly, between 7% (C scheme) and 
19% (J scheme). The i^-factors, for fi = yfs, are also very different, from one algorithm to 
another, again with the minimum corresponding to the C scheme (K = 1.07) and the maximum 
to the J one (K = 1.20). All these values are however lower limits. In fact, as fi is decreased a s 
increases, hence the relative size of the NLO effects grows larger too, in each case. 

As an estimator for the W mass we use the 'average' mass, M ave , defined as follows. Out of 
the three possible combinations of pairs of jet-jet systems, we choose the one for which the two 
reconstructed W ± masses, M Rl and M R2 , minimise 

AM = \M Rl - M w ± | + \M R , 2 - M w ± | (6) 

and then define 

M ave = ^(M Rl +M R2 ). (7) 

This variable has been extensively used since, at tree level, the difference between M ave and the 
average between the two W ± masses that one would reconstruct if the quarks could always be 



7 Note that by restraining /i to values higher than the hadronisation scale Qo, which is of order 1 GeV, a 
perturbative analysis is in principle always justified. However, too low a value of /i would imply a very large 
a s , in turn rendering the fixed order predictions unreliable. As a compromise, we will be considering /i- values 
well above Qq in the reminder of our study (say, 35 GeV and above). For scale choices in the interval 35 GeV 
< /i < 175 GeV, the strong coupling constant varies over the following ranges: 0.162(0.134) < a s < 0.123(0.105) 
at one-(two-)loop level. 
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Figure 2: The zero-th (LO) and first order (NLO-LO) a s components of the NLO differential 
distribution in the 'average' mass (as defined in the text) for e + e~ — > W + W~ — > udsc (CC03), 
the latter with QCD scale /i = y/s, as given by four different jet clustering algorithms, at 
\fs = 175 GeV. (The values of y cut are chosen such that the LO rates are approximately equal 
for all schemes.) Notice that the dashed and dotted blue-lines are visually indistinguishable. 
Total hadronic rates are obtained by multiplying those above times four. 



paired correctly is very smallQ. 

The NLO differential spectra in M ave show a shift towards low mass values with respect 
to the LO case which depends upon the jet algorithm being used and its y cut value. (The 



generation of this low mass tail at NLO was already observed and discussed in Ref. |p3| , where 
its consequences for a determination of M^± were described in details.) 
The NLO distribution in a generic mass M is made up by two terms: 

da NLO /dM = A(M) + a s B(M), (8) 

one proportional to the zero-th power of a s (denoted by A: the Born term) and another to 
the first power (denoted by B: the first order correction). They are plotted separately in 
Fig. H] for the mass difference M = M ave — M w ±, labelled as 'LO' and 'NLO-LO', respectively. 
Numerical values of the two functions for some selected 'average' masses are given in Tab. [I] 
(where M = M ave ), both normalised here to the total Born cross section (i.e., the integral over 
M ave of the LO curves in Fig. |2|). The scale of the strong coupling constant is still set to \i = y/s. 

8 One could consider more sophisticated approaches, but this is beyond the scope of this paper. 
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Figure 3: NLO differential distribution in the 'average' mass (as denned in the text) for e + e _ — > 
W + W~ — ► udsc (CC03), for three choices of the QCD scale fi, as given by four different jet 
clustering algorithms, at \fs = 175 GeV. (The values of y cut are chosen such that the LO rates 
are approximately equal for all schemes.) Total hadronic rates are obtained by multiplying those 
above times four. 



However, we are concerned here with the fact that the actual shape of the M ave — Mw± 
distribution at NLO depends upon /j (this was set to yfs as default in f23j), through the choice 
of both the jet finder and its resolution parameter. Ultimately then, so will do the value of 
Myy± extracted from that distribution. To study this effect, we plot in Fig. |3] the differential 
distribution of the quantity defined in eqs. (§D-([7p, now for three different choices of /i, e.g., 
35, 100 GeV and y/s, for our default choice of jet algorithms and y cu t's. There exists a visible 
variation with /x; besides, the previously observed dependence on the choice of the algorithm 
and y cut persists at different /j's. 
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Table 1: The A(M ave ) and a s B{M scvc ) components of the NLO cross section normalised to the 
Born rate, see eq. (||), for representative values of M ave , as obtained by our default jet clustering 
algorithms and resolutions, with a s = 0.123 evaluated through one-loop order at the scale 
/i = a/s = 175 GeV. Recall that the input value for the W ± mass is 80.430 GeV. 
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Figure 4: The NLO mass difference 5M =< M ave — M w ± > as a function of a s calculated at 
one-loop in e + e~ — > W + W~ — > udsc (CC03), for three choices of the QCD scale /x, as given by 
four different jet clustering algorithms, at y/s = 175 GeV. 



In order to quantify the impact of the /i-dependence of the W + W~ — ► 4 jet rates on the 
W ± mass, we have collected in Tab. ^| the values obtained for the mean deviation from M w ± 
of the reconstructed 'average' mass M ave , as predicted by our default jet clustering algorithms 
and separations. Both at LO and NLO, the difference < M ave — M^± > is negative, as already 



observed in p3fl . Besides, by comparing the higher order predictions obtained with \i varying 
from 35 to 175 GeV, one may notice that systematic uncertainties on M w ± could turn out be 
very large in the end, since < M ave — M w ± > can be as large as 500 MeV (in the L scheme). 
The same data are plotted as a function of a s in Fig. |4j For all schemes the a s dependence is 
essentially linear and the shift for different ranges in a s is readily evaluated. 

In order to estimate more realistically the systematic uncertainty on the determination of 
the W ± mass induced by the unknown scale fi, we perform a MINUIT |36| fit on the LO and 
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< M avc - M w ± > (GeV) 


Algorithm 


/i = ys/o = oo bev 


/l = i/S = 1UU Lrev 


, , /T. 1 *7K (~~* r^\T 
fi = a/S =1/0 VjreV 


J 


-0.34 


-0.34 


-0.34 




-1.91 


-1.67 


-1.57 


D 


-0.28 


-0.28 


-0.28 




-1.71 


-1.49 


-1.40 


C 


-0.28 


-0.28 


-0.28 




-1.94 


-1.66 


-1.55 


L 


-0.35 


-0.35 


-0.35 




-2.49 


-2.13 


-1.99 



Table 2: Mean difference between M ave and My/± as obtained from (some of) the spectra in 
Figs. |^-|3|. First line is for LO results (these do not depend upon fi), second line is for the NLO 
ones. Recall that the input value for the W mass is 80.430 GeV. 



NLO distributions, with a fitting function of the form 



f{ m )= C lT~2 1\2 i 2 2 +#( m ) ( 9 ) 

(m 2 - ci) 2 + 44 



where the term g(m) is meant to simulate a smooth background due to mis-assigned jets induced 
by the clustering algorithm. For the latter, we adopt two possible choices 

( c 4 + c 5 (to - c 2 ) + c 6 (to - c 2 ) 2 , 
^(to) = c t (10) 

I C 4 l+exp((m— C5)/ce) ' 

that is, a three-term polynomial and a smeared step function (motivated by the kinematical- 
limit shoulder at large masses and on the same footing as in Ref. [|l3j). Notice that in eq. @ 
we have implicitly assumed a Breit-Wigner shape characterised by a peak height c±, a position 
c 2 and a width c 3 , corresponding to the normalisation M w ± and Tjyi, respectively, of the 
distributions in Figs. To first approximation, the difference between the values of the 

coefficient c 2 as obtained from fitting the above curves is then a measure of the typical size of 
the systematic error that we are investigating. Obviously, more sophisticated fitting procedures 
could be adopted, possibly yielding different results for M w ±. However, it should be clear from 
a close inspection of the plots in Fig. and the estimates in Tab. [| that varying /i over any 
reasonable interval would result in mass shifts comparable to or larger than the uncertainty 
reported in (|l|). 

Tab. [3] reproduces the results of one of our fits. Whereas the values of some of the parameters 
(such as the height h and the width T w ± of the distributions) depend sensibly on the choice of 
the mass interval used for the fit and/or the form of the background, the values extracted for 

9 This is related to the input normalisation for MINUIT and is irrelevant to our purposes. 
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^ = = 35 GeV 


V = V 


>s = 175 GeV 






M w ± 


r w ± 


h 


M w ± 


T w ± 






(GeV) 


(GeV) 




(GeV) 


(GeV) 




Pol 


ynomial background 


T 

J 


517.757 


80.464 


3.237 


575.412 


80.454 


2.828 


n 


480.344 


80.298 


3.163 


539.973 


80.343 


2.788 


c 


422.196 


80.211 


3.450 


486.111 


80.292 


2.930 


L 


453.062 


80.085 


3.801 


529.413 


80.217 


3.142 




Smeared step function background 


J 


516.169 


80.451 


3.213 


601.842 


80.472 


3.008 


D 


521.475 


80.286 


3.450 


549.312 


80.332 


2.890 


C 


454.114 


80.226 


3.619 


517.671 


80.320 


3.163 


L 


456.933 


80.154 


3.791 


515.275 


80.221 


3.127 



Table 3: Fits to (some of) the M ave spectra at NLO in Figs. First three columns are for 
\i = \/s/5, last three are for \i = y/s, with y/s = 175 GeV. A Breit-Wigner shape is always 
assumed, supplemented by a three-term polynomial (upper section) or a smeared step function 
(lower section) to emulate the intrinsic background. Here, we have fitted the M ave distributions 
over the mass interval 75 to 85 GeV. Recall that the input value for the mass is 80.430 
GeV. 



Myy± at NLO are remarkably more stable^. The parameters in the table should be taken as 
representative of the qualitative features of all fits that we performed. From there, one notices a 
strong dependence of the fitted M w ± value upon the jet scheme (for a 'fixed' /i), with variations 
of up to almost 380 MeV (between the J and L schemes, when /x = 35 GeV and assuming a 
polynomial background). 

However, all effects discussed above are well understood, since they are merely of kinemat- 
ical origin (the different handling of gluon radiation by the various jet clustering algorithms). 
As already stressed repeatedly, it is the systematics associated with the choice of \x that is be- 
yond theoretical control. It turns out that such an intrinsic uncertainty of the fixed-order QCD 
predictions can be rather large in the determination of the mass: compare the values for 
M\y± as obtained in the NLO fits and given on the left of Tab. |3| to those on the right. The 
differences between the reconstructed M\y± values for \i — 35 and 175 GeV (which, hereafter, 
we denote by 5nlo^w±) can become as large as 130 MeV (for the L scheme, in presence of a 
polynomial background). The J scheme seems to be here the least sensitive to \i- variations (a 
10 to 20 MeV effect), with the D and C ones falling in between. Notice that, in this exercise, we 
have restrained ourselves to values of /i not smaller than \/s/5, precisely in order to avoid the 
mentioned logarithmic effects induced by a choice of fi too close to the jet-scheme-dependent 
emission scale y/y cu ts (see also Footnote 7). Indeed, for /i in the above interval, we have found 

10 For reason of space, we do not reproduce here the values of the coefficients C4, c 5 and Cq which characterise 
the background. 
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that to change the values of y cut (still maintaining the typical four-jet separations used in ex- 
perimental analyses) does not affect our conclusions. In fact, the latter do not change either, if 
one adopts a different recombination procedure of the cluster momenta. 



3 Conclusions 

In summary, we have verified that theoretical errors arising from the use of different jet clus- 
tering algorithms in treating the fixed-order 0(a s ) corrections from perturbative QCD, in the 
prediction of experimental observables used for the extraction of the W ± boson mass in the 
hadronic channel e + e~ — > W + W~ — > 4 jets (CC03) at LEP2, can be competitive with similar 
systematic effects that could be induced by non-perturbative dynamics, such as CR and BECs, 
e.g., as predicted in the hadronisation model of Ref. f37| . In particular, using various jet defini- 
tions, and spanning the scale of a s between My/± and 2M^± approximately, we have shown that 
these uncertainties on My/± can be of order 100 MeV, hence larger than current experimental 
assumptions on the size of the theoretical error. 

We also have verified that the inclusion of irreducible background via CC11 diagrams has lit- 
tle impact on our main conclusions, so has the incorporation of ISR effects. Similarly, a different 
treatment of the widths in the resonant propagators (the fixed-width scheme was adopted here) 
has negligible consequences for both CC03 and CC11. A different choice of \fs (at fixed 2/ C ut's, 
or vice versa) yields similar estimates of 5^oM w ± to those given here. Also, if one enforces 
typical W ± mass reconstruction cuts, say, \M Ri — M w ± \ < 5, for i — 1,2 and 10 GeV < 5 < 30 
GeV, see eqs. (§)-(0)) typical values of 5nlo-^w± remain in the above range, despite the effects 
on the actual event rates can be dramatic |23| . 

Finally, effects due to the kinematic interplay between jet clustering algorithms and PS 
(including hadronisation) were not in the original intentions of this study, as they have already 
been addressed in Ref. ||13|| . Whereas the latter can be estimated in the context of an event 
level MC analysis, those considered here are intrinsic uncertainties of the theory. The results of 
our present analysis point to the fact that such perturbative QCD effects may not yet be under 
control, at least in the context of mass determinations from the hadronic data samples 



collected at LEP2. Taking also into account the results of Ref. |T3[ in the same context, wherein 
the systematic uncertainties in the reconstructed M\y± value due to the dynamics involved 
beyond the hard scattering processes (as obtained in HERWIG and PYTHIA) were often found 
to be somewhat smaller than 100 MeV ( more in line then with the experimental estimates 
discussed in the Introduction), one may conclude that a thorough reassessment of the theoretical 
systematics entering the e + e _ — > W + W~ — > hadrons channel is in order, given the importance 
that the precise knowledge of the mass has in constraining the properties of yet undiscovered 
particles, such as the Higgs boson mass. This will require the availability of a MC event generator 
based on the NLO matrix elements (both real and virtual) used in this analysis (properly 
interfaced to the subsequent PS and hadronisation stages), which is under construction in the 
HERWIG environment [P5| . It will eventually be the implementation of the more sophisticated 
selection methods used by the LEP collaborations (as opposed to the simpler ones illustrated 
here) in the context of a such a NLO-MC event generator that will finally assess the uncertainty 
range on M w ± still allowed by the most up-to-date theoretical and experimental instruments. 
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